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The General Circulation of the Atmosphere. 

Bt Francis R. Sharpe. 



§ 1. Introduction. 

The difference in temperature between the equator and the poles of the 
earth tends to produce a movement of the atmosphere from the poles towards 
the equator at the earth's surface, and thus set up a steady circulation. The 
problem has been investigated by Ferrel* and Oberbeck.f The former did not 
attempt to make a complete analytical solution. The latter assumed that the 
atmosphere was incompressible but that it satisfied Boyle's law for a perfect gas, 
and he surrounded the atmosphere with a spherical boundary of undetermined 
height. In the present paper the known hydrodynamical equations for a viscous 
gas are first determined, using polar coordinates, by considering the flow of 
matter and momentum through a small volume. By finding the flow of energy 
another equation is obtained, which is used instead of the ordinary adiabatic 
assumption. This equation is equivalent to the equation of energy in the 
kinetic theory of gases but is here obtained independently of that theory. 
Making use of the fact that the height of the atmosphere is small compared with 
the earth's radius, an approximation to the motion is found for a non-rotating 
earth. The efiiect of the earth's rotation on the east and west motion is next 
considered, and then the resulting modification of the pressure gradient. Finally 
the solution which is obtained is compared with the observed facts. 

§ 2. Coordinates and Variables. 

In discussing the general circulation of the atmosphere on a spherical earth, 
it is convenient to use polar coordinates; namely: r, the distance of any point 
of the atmosphere from the earth's center ; 6, the colatitude ; and ^, the longi- 
tude. The variables which we wish to determine as functions of r, 6, 4>, and t, 
the time, are : p, p, r, the pressure, density, and absolute temperature ; and 
u, V, w, the components of the velocity. 

* Ferrel : Recent Advances in Meteorology. 

•)• Oberbeck : Papers on the Mechanics of the Earth's Atmosphere, translated by Cleveland Abb6. 
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The atmosphere is supposed to be a perfect gas, so that 

p = Rr^, (1) 

where i2 is a constant. Five other equations are necessary to determine the six 
variables jp, p, t, u, v, w. They are derived in the next three sections by con- 
sidering the flow of matter, of momentum, and of energy. 

§ 3. The Flow of Matter. 

Consider the flow of matter through the six faces of a small volume 
dr . rd$ . r sin dd^ bounded by the spheres r, r + dr, the cones 6, 6 ■{■ d6, and 
the planes ^, 4> + d^. The rate of increase of the matter in this volume is 

-^{dr . rd$ . r sin $ d^) = — -^ (pu . rdd . r sin 6 d^) dr 

— ofl (pv • *■ sin B d^ • d/r) dd — ^ {pw . dr . rdd) d^. 
Dividing by dr.nHS.r ain$d^, we have the continuity equation: 

^ + ^3-(''^P") + ^^-o^(p«'sin0) + — ^-^o-(pM') = 0. (2) 

dt r^dr^ '^ ^ rsmOdO^^ ^ rsmdo^^'^ -^ ^ ' 

§ 4. The Flow of Momentum. 

Consider the flow of momentum through the same small volume. Momen- 
tum being a vector quantity it is necessary to take account of the changes in the 
direction of the components of the velocity as we consider the different faces of 
the polar element of volume. Thus u turns through the angle dd on the plane ^ 
and through sindd^ on the cone 6, while v and v + dv meet at the angle 
cosdd^ on the earth's axis. 

The normal stresses on the faces P, Q, B, the tangential stresses S, T, U, 
and the external force F also generate momentum; moreover the stresses 
change their directions as we consider different faces. The rate of increase of 
momentum in the direction of r is therefore 

-^ (pw) dr . rdd .rsindd^=^ 

— ^\{pu.u + P)r sindd^ . rdO} dr —^\(pv.u-{-U)r sin dd^ . dr\ dd 

_ J_|(ptp.M+ T)rdddr\d^ + {pv.v -\-Q)r sin Qd^dr .dQ 
+ {pw .w ■\- B)r dB d/r sinB d^ -^ p F^dr .r dB . r sinQ d^. 
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Hence we have the equation of motion : 



r^dr 



r sin % dd 



+ 



r sin 6 d^ 



(p«,.^+y)_P"-^ + g + P^-^ + ^_pj;==0. (3) 



Similarly, 



dt 



and 



Ir(p^) + ;:A;.Kp«-^ + f^)'-n + r^T^.-Kp^'-^ + e) 8in0f 



r^dr 



+ 



rsin 



^(pi«.« + ^)- 



r sin $ dd 
pv.'U+U—{pw.w + Ii)cotd pj?.^ — Q (4) 



■p^«: 



^(P^') + A: Kp« . w' + ^) '•n + ::r-|^ Upt' • ^^ + '^) sin e^ 



r^9r 



r sin $ dd 



L 3 /' j^ r>\ 1 P«'.'M + 7'+(pt».tJ + ^)cot0 ij. _n /'K\ 
<- ^.;r,fl;^^ (p^-^ + -^) + - 1;. P^* = 0. (6) 



rsin 034) 



9p 



If p^ is eliminated by means of (2), these equations can easily be identified 
with the ordinary form of the equations of motion in polar coordinates. 

In Stokes' theory of viscosity the stresses P, Q, B, S, T, U are linear 
functions of the strains e,f, g, a,b,c; namely: 

P = p + %fiS-2iie, 

B=p + %(i8— 2ftg, 
S = — [I a, 
T=-(ib, 
U= — [lo; 

where the dilatation 8 = e+/+g; the three extensions are 

du 



(6) 



e ^ 



dr 



, dv u^ 

•^ " rde "^ r ' 

dw u 

^ ~~ r sindd^ r ' 



(7) 
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and the three shears: 



dv , dw w , A 

« = — • nc,^ + -^a ^^^ ^> 

r sm vo^ rou r 

7 du , dw to 

r sin 6d^ dr r ' 

du , dv v_ 

rdd dr r ' 



(8) 



On substituting the values of the stresses in the equations of motion we 
could obtain them in the form which is generally used, but it is more convenient 
for our purpose to leave them in their present form. 

§ 5. The Flow of Energy. 

By considering the flow of matter and momentum we have obtained the 
equations of continuity and motion. In an exactly similar way the kinetic 
energy per unit of volume is \^{v? -\-v^ •\- vF), and like p it is a scalar quantity. 
The energy in the form of heat is (7pT, where G is the capacity of the atmos- 
phere for heat. The rate of conduction of heat is — x grad t, where k is the 
conductivity. The stresses P, §, iB, 8, T, U, and the external force F also 
do work and thus contribute energy. Taking account of all these different forms 
of energy, we find for the rate of increase of energy in a small polar element 
of volume : 

i- lip («*^ + »* + ^) + C'p'^l <^ -rdS.r Bindd^ = 

— ^[ hp{u^ + v' + v^)u+ Cpru + Pu + Uv + Tw — xplrdSr amOdp^dr 
_ ^rUp(tt2 + ^2 + «,8)tj + Cptv + Uu + Qv + Sw — x^ir smed^dr~\de 

— ^\hn(u^+ir' + w^)w+Cp'tw+Tu + Sv + Bw — x—^^^^lrdedr~\dA 
3^ LI rsmdd<p) J 

+ p{FrU + FeV + F^w)drr dOr sin Bd^. 

Using (2) to eliminate ^ , substituting for the three components of F from (3), 
(4) and (5), and using (7) and (8), we find 
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Gf^-itPe+ Qf-^Bg + Sa+Th+Uc 

9 /^,,^«9t\ a (:c^ nine) ^ ( :c ^ \r=^. 

r^dr\ dry rsin6dd\ rd$ / r8indd(p\ raindd^/ 

On substituting the values of the stresses from (6) this equation becomes 

r^dr\ drJ rsin039\ rd% / r%va.Qd^\ rsmdd<py 

This result agrees with the equation of energy in the kinetic theory of gases, 
but is here derived independently. 

Replacing p and S by the use of (1), (2) and (7), we obtain the more 
convenient form 

7^dr\ or/ ranidad\ rod v ra\TxQo^\ rsm09^/ 

The conductivity and viscosity of the atmosphere are both small ; hence, 
if the motion is so slow that we may neglect squares of velocities, (9) becoraes 

^ dt dt 



Integrating, we find that 
which gives, from (1), 



<r-cc PS 



p cc p o 



R + C 
t 



the adiabatic law. 

Since H is about . 4 G, we have 

p=zmp^*, (10) 

where m is a constant. 

§ 6. Adiabatic Equilibrium. 

When the temperature is a function of r only, the equations of motion 
become 

^ + .P = o, 
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or, from (10), 



1 . 4 wp* J- + gfp = 0. 



Integrating, we have 

1.4 
.4 

1.4 



that is, 



mp* -\- gr = constant ; 
Er -\- gr = constant. 



.4 

At the surface of the earth r is a and r is Tq ; therefore 

1.4 



.4 
Hence t is when 



B{t — %) + g{r — a) = 0. 



r = a + ^Bto. (11) 

Substituting numerical values, we find that the height of the atmosphere is 

a, which will be denoted by ah, where h = — — . 

220 ' 220 

If we now put r = a (1 + ^ — as), so that ax is the distance from the top 
of the atmosphere, then r is a linear function of x which vanishes when a; = 
and has the value -To when a; = A. 

Hence, we have 

^ =<%)■■ ('^) 

and therefore, from (1) and (10), 



P 



=<U <>" 



and 



p = ^(f)'. (H) 

A relation between the constants is furnished by (1) and (11); namely, 

a = -^^. (15) 

2 a Apo 

8 
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§ 7. Circulation on a Non- Rotating Earth. 

Assume that the temperature of the atmosphere decreases symmetrically 
from the equator to the poles and diminishes uniformly from the surface upwards. 
Then, instead of (11), we have 

T = -^(<ro — r^cos^e), (16) 

n 

where -^ is small. 
■To 

To satisfy (1), and from the form of the equations of motion, we therefore 

take instead of (13) and (14) : 

i'=i'o(|) + PzGOs^B, 

h 
where 

Pi _ pa _ ^ 

\^ To 



p = Po(|) + pgcos'^e, 



/xY /xY 

^<h) ^<h) 



(17) 



(18) 



Equilibrium is no longer possible, because, from (6) and (17), 

and (4) is therefore not satisfied. 

§8. Solution of the Continuity Equation. 

When the motion is steady and symmetrical, the continuity equation takes 
the form 

^ (P^") + -^A^«(P^«^'^^) = <^' 



i^dr r sin B dd 

or approximately, from (14), 

/x\i { du , d , -ml 6 /aj\i „ 

( — ) \ — ^ — h -. — TT^ (v sm 0) y r ( -^ ) w = 0. 

\h) \ ax sin^ae^ ') 1h\h) 



Hence 



9!i + AJi = -_^ (r;sin0). 

8a; ^ 2 » sin^ae ^ ' 
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The solution of this equation, which corresponds to (17), is 

« = Cr(i-3cos«0), 1 

« = Fsin0cos0, J ^ ^ 

where U and V are functions of x which satisfy 

§ 9. Approximate Form of the Equations of Motion. 

Since a; is a small ratio which lies between and h (that is, between and 
^^), we see from (20) that U is small compared with V. Therefore, if we 
retain in the equations of motion only the most important terms, we find 

^+9f = 0, (21) 

From (15), (17) and (19) these equations become 

^ = LmL (23) 

dx 2 hpo 

and 

Since the viscosity varies as the .76 power of the absolute temperature, the latter 
equation may be written in the form 

— i(-(l)'II)- (-) 

§ 10. Solution of the Equations of Motion. 

Finding first jp from (18) and (21), next Ffrom (24), and then ZJfrom (20), 

we have 

U = Ax-i + BxV^ + GxV- log x + Dxi + Ex, ( 25) 

where A, B, D and E are arbitrary constants, but C must be so chosen that (18) 
is satisfied.- 
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§11. The Boundary Conditions. 

At the top of the atmosphere (where x is zero) it is assumed that u is zero 

and that the tangential stress /K^r— is zero; hence A and D are zero. 

or 

At the surface of the earth u and v are both zero ; hence 

U=bQ^ jbV log ^ — ar^t^ + AV cc) . (26) 

Substituting this value of Z7 in (18), we find 

^^^^960^0^3 /^S^y iog^_«,¥H- 7t'^a:Vl _ 3 cos«0), (27) 

107217AV^o'z^o^ 4 ^h ^ J^ '' ^ ' 

896a/)oT2 /627 ;. , a; 7 ,, , 7 ,jl.\ . /, ^ ,„„x 

w = i^r-^ — ( ^^ log -r a; ^ + -h^' ) sm cos ©, (28) 

107217pV^oT'oV 16 ^A 2 ^2 / ' ^ ' 

p=p,(fy{l + (Ilogf+?^)eo«'<.}. (30) 

§ 12. East and West Motion on a Rotating Earth. 

When the rotation of the earth is taken into account, it is convenient to 
denote by w the east and west velocity relative to the earth so that the actual 
velocity is w + er sin 6, s being the earth's angular velocity about its axis. 
Substituting this value in (5) and retaining only the most important terms, 
we have 

= - 2epu sine- 2^pt, cosO + -|-(^^) ; 
hence, using (19), 

|-Ca;i|^) = ^^^°^ £^8ine(l - 3 cos^e) + Fsin^ cos'^eh 

The term containing U is small, compared with the term containing F, except 
near the equator, where cos 6 is small. Substituting from (27) and (28) for U 
and V, integrating twice and remembering that at the top of the atmosphere 
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dw • 



(where x is zero) the tangential stress /* 3- is zero and that at the surface of the 
earth (where a; = A) «> is zero, we find 



w 



1072l7M^roLl37 ^h 19. 37=^ ^171 

-^^zf^'^l^^H^-^^^os'^) [(32) 

xV + A AJ^a;¥ - -^ h^] sin e cos^ ^l . 
15 4335 J J 



134 ^A 289 



§ 13. Pressure Gradient Due to the East and West Motion. 

When the rotation of the earth is taken into account, the equations of 
motion (21) and (22) become 



^ dp . . 

= — ^ — 9P + iepw sm0, 
or 

dp 



dv\ 



' = -^ + '^P"*'''^ + ^(''' drr 
We must therefore replace (17) by 

p = po (|-) + Pa cos^ + p^ cos* e, 
P=Po(^j) + Pi cobH + Pi cos* e, 

in which, from the law of gases, as in (18), 

Pi _ P« 



(33) 
(34) 



(35) 






(36) 



Again, (19) must be altered to 

M=: U{l — 3cos^e)-\- Ui(SeosH — 5 cos* 6), 
V =z V sinB cos B -\- Vi sin 6 cos^ 6, 

where (from the continuity equation) 



dx 2 X 



V,. 



(37) 
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From (34) we see that p^ is of the same order as p«c in x, while V^ is of 
order two higher in as; and therefore, from (87), C4 is of order three higher in x. 
Hence (33) becomes 

= — -— — go. 
dr ^^ 

Using (36) and (15) we therefore find, on integration, 

Pi — P^'^, 
where P is an arbitrary constant which will be subsequently determined so that 
the boundary conditions are satisfied. 
Substituting in (34), we have 

0= APxiBine<iOs^e+ 2aepw cos0+ ^ -^ ^a;||^^ sin cos^e. 

Substituting the value of w from (32) and retaining only the second term 
in that equation, we have 

_ , 3684a*e^2'oPo'^3 f 19 in x 53 u , 4 ,19 ^, 361 ,1, 5) 
107217 AV§ -To (34 ^h 289 ^15 4335 ) 



am ax \ ax / 



Integrating with the boundary condition that, when x is zero, the tangential 
stress /w-o^ is zero, we find 

8 „ , 3584a*62j9oPo'^3 f 19 la , x 1595 13 , 16 7,19^29 

= Pxi ^"' — ? i X loff -5 X + ft 5 x~i~ 

Q 107217 AVi-^o 1 408^ "^ A 288.289 ^435 

722 iT.iJ , «o 3 3F4 
30345 S ahi dx 

Dividing by xi and integrating with the boundary condition that, when x 
is h, Vi is zero, we have 

n_32p,i, ,,9>, 3584a^g^yop§ra < 19 „„«. 1 _ „ «! 93559 y 

^ -rn^i^^-^'^)- io72177iV^To' 14998''^ ^^^ A - 72. 17^49^ ''^ 

, 32 X 14 15. 2888 1.11 15 87, 822, 265 T4_9f • Wo rr 

+ 4525^"^^"-465175^""^- 52200. 17^49^ ^^"1 +ili ^^- 

Substituting the value of V^ in (37), multiplying by a;t, integrating with the 
boundary condition that CT^ is zero when x is zero, and, finally, putting x equal 
to h, we have, as the condition that Ui is zero, the equation to determine P : 
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171 V33 7/ 107217 AV§ To i 4998 63** 72.17^49^63 



("51^0 

87^ 

6525 11 455175 29 52200. 17^49^ 



, 32 J^ _ 2888 ^ 87,8 22,265 21 



Hence 



107217 /«§To ' 



and therefore 

/a;\H, , T2/7, „x , 119\ „o , 140.7a*e^»oPo''^2^^ d/^) /^ n 

P=P<-h) 1^+^(2^"^^+ 9T)«^^^ + --mir^f^ ('') 

§ 14. Discussion of the Results Obtained. 
The vertical velocity u^ has a maximum value 

'i^iapQt^h^ 

107217 jUo'^o 
when 05 is about .55 A. 

The southerly velocity v^ has a maximum value 

1023ap„T2A^ 
107217 f/on 
when a; is about .886 A; and it changes to northerly when x is .73 A. 

It is seen from (32) that the velocity is easterly except near the equator, 
where, on account of the stnallness of cos 0, the second term is important and 
the velocity is westerly. 

The pressure p^ has its maximum positive value, 

119j?oT2 

99 To ' 

when X is h. It vanishes when x is .71 A and has its greatest negative value 
when X is .53 A. 

§ 15. The Observed Facts in the Atmosphere. 

It is known that* the temperature at the surface of the earth is 

8.5 — 20.95 cos^e, 
or, in our notation, 

* Ferrel : Recent Advances in Meteorology, p. 462. 
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r = 303 — 42 cos^e, 
= To — Tg COS* 6; 

so that — * is not very small. It is also known that * 

^ = 758 + 31 cos^e — 61 cos*0. 
Putting 33 = A in (38), the value obtained was 

^-P' + 99 To """^ ^ 107217 ^gTo ^^^ ^' 

hence our value for p^ appears to be too large, although its value decreases 
very rapidly as we ascend vertically. As regards jj^, if we substitute numerical 
values we find too large a value. The explanation is probably that, on account 
of the large amount of eddying motion in the atmosphere, the ordinary ex- 
perimental value for ^ is smaller than its effective value in the atmosphere. 
Moreover the motion we are considering is unstable except for small values of 

^^ and large values of /w. There appears then to be only a qualitative and not 

a quantitative confirmation of the solution obtained. 

As regards the velocities, it is agreed f that the velocity at the equator is 
westerly, and in the middle latitude easterly, and increasing with the altitude ; 
but there is no evidence of a northerly motion in the upper atmosphere as theory 
would lead us to expect. 

CoBNELi. Univbbsitt, July, 1908. 

* Overbeok : loc. cit. 

t Hildebrandson : Brit. Assoc. Report, 1908. 



